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Abstract 

We employ QCD sum rules, implemented with two numerical algorithms already 
tested in two different channels of Charmonium, in order to predict masses of resonances 
just above the ground states in the 3 P channels. We find that such masses are above 
the threshold of open charm. We calculate also the partial decay widths of the ground 
states into light hadrons and, for even spins, into two gammas; we find consistency 
with data. 
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A lot of experimental and theoretical efforts have been performed in order to study 
hadrons formed by one or more heavy quarks. In particular, as regards spectroscopy predic- 
tions, essentially three different approaches have been adopted in the literature: i) potential 
models]!], [|; ii) lattice calculations [|J; in) QCD sum rules[|5], |], [7], Concerning heavy 



quarkonia, a systematic study of spectra has been done on Bottomonium|| [10| , but not on 
Charmonium. We think the gap can be partially filled in by means of power moment QCD 
sum rules f|, || || [TT], |j~2f , which, incidentally, in the case of Charmonium appear more suit- 
able than exponential sum rules|7|, §, [H], [L3|]. These have been already used successfully, with 
the help of two rather complex numerical algorithms, in the 1 Sq and 1 P\ channels |TT|, |i~2|l of 
Charmonium, in order to determine the mass of the first excited state - be it a bound state 
(rj' c , 1 S channel) or a resonant state (h' c , 1 Pi channel) - and, as a byproduct, some partial 
decay widths of ground states. In this letter we extend our investigation to the 3 P-channels 
of Charmonium; as a result we find that they have no bound states besides the ground one. 
Moreover we calculate the partial decay widths of the ground states into light hadrons and 
(for even spin resonances) into gammas, finding consistency with data. 

We recall shortly the power moment QCD sum rules, together with the two numerical 



algorithms, which we have elaborated in preceding papers|TT], |12| . We assume a dispersion 
relation for the Fourier transform of the two-point function of the current (or density) of 
the charmed quarks with the quantum numbers of the channel considered. We consider 
spacelike overall momenta, q 2 , therefore we may develop the two-point function according to 
a Wilson expansion - a nontrivial one, including the contribution of the lowest-dimensional 
gluon condensate. Then for a given channel T and for any positive integer n we get [|TTH 



^[1 + a s al(0 + <K£(0] = \\ ^Q^ ds, (1) 



where Q 2 = —q 2 , a s = a s [A{m° c ) 2 + Q 2 } is the running coupling constant (a s [4(m°) 2 ] 
Q 2 

4(m») 2 ' 



a s ~ 0.3), £ = tAw, m" = m c (p 2 = -m 2 ) 



*=36(i <0| ?W 0> < 2) 
and G a is the QCD strength tensor field. The l.h.s. of eq. ([1]) coincides, up to an n- 
dependent factor, with the n-th moment of the Wilson expansion. Concerning the spectral 
function, which appears at the r.h.s. of eq. ([TJ), for heavy quarkonia it may be parametrized 
as follows: 

Imlf{s) = ^ £ ^ « -*) + ^(1 + f M« - so), (3) 

where m ; is the mass of the i-th resonance of the channel and g>j the coupling constant of 
that resonance to the current (or density) relative to the channel, while sq is the threshold 
energy squared of continuum; lastly <Tq is a positive integer, which we set equal to 1 for 
P-channels[Q. Owing to duality, we are free to assume N = 1 or 2 in the sum that appears 
at the l.h.s. of eq. (|3|), including the other resonances in the continuum contribution. If we 
take N = 1, from eqs. ([!]) and (|3|) we deduce the mass of the ground state, i. e., 

mi = ^~Q 2 , Mf = Mj- ,i ( ' + ^ n , (4) 
M n+1 8n 2 n(s + Q 2 ) n 

where M v n = A$[l + «X(£) + $ i s the l.h.s. of eq. ([![). If, instead, we take N = 

2, we get the mass of the first excited state, i. e., 

Equations (f|) and @ yield respectively mi and m2 as functions of n, Q 2 , a s , $, and 
So, the last parameter varying from channel to channel. m 2 depends also on g x , which is 
not known for most channels. Preceding analyses yield $ = 1.35 1CT 3 || [/J, a s = 0.3@, || 
and mP c = 1.268 GeV^Ill]]. Moreover, as we have already shown [11|, we can fix criteria for 



determining the other parameters. 

In principle the value of sq that appears in formula (^) could be different from the one 
used in formula (|5]), since in the latter case we exclude the contribution of irt2 from the 
continuum of the spectral function. However eq. ([3]) is a simplified parametrization of the 



real spectral function, which presents narrow resonances at low energies and therefore has 
a greater weight at low s in the integral at the r.h.s. of eq. (|]). Then s results to be a 
decreasing function of n and Q 2 . We can set Sq ^ — Sq \ provided n\ ^ 112 and Q 2 ^ ^ Q(2)> 
the index 1 (2) referring to the option N — 1 (2). This self-consistency condition is always 



respected, both in the present analysis (see below) and in preceding ones [|TT| , |12 , [14], |15|1 . 

In order to fix n and Q 2 , first of all we consider the dependence of m^j) on n at a fixed 
Q 2 - rrii(n) present minima at given TTj, % = 1,2, and therefore plateaux around such values 
(see e. g. fig. 1). Furthermore the amplitude of each plateau, which can be measured by 
the inverse of the quantity 



Di = mi{m - h) + mi(ni + k { ) - 2rrii(ni), (6) 

(ki = 2, &2 = 1), varies with Q 2 . We choose so as to coincide with locations of minima 
of Di (see figs. 2). In particular we pick up for £1 the lowest local minimum (at least for not 
too large £) of Di, provided the corresponding mi is stable with respect to the parameters m° c 
and Sq,; as we shall see below, the stability condition in the 3 Pi channel requires a particular 
care. As to £2, we pick up the lowest local minimum of D2 which does not coincide with any 
minima of D\. 

We determine s$ by equalling mi to the mass of the ground state, which is generally 
known experimentally. In order to determine gi, we proceed as in the vector channel flT[ 



considering the graph of D2 (that is, of the value of D2 obtained by selecting the minimum 
in the way just described) versus g\. The value corresponding to an oblique inflexion (fig. 
3) is assumed to be the correct one, provided it is stable with respect to small variations of 
the parameters. This in turn allows to determine m2 by the formulae exposed above. 



Table 1: Results of the analysis in different channels of Charmonium 





m 2 (MeV) 


9i 


so(GeV) 2 


n 


Q 2 (GeV) 2 


"Pi 
3 P 2 


4097 + ^ 
4327l| 
4466±H 


1 q Qo+u.34 

lO.OO_ .07 

16.0+g-f 

ii.82t°':| 


90 7 +U.12 
zo - ' -0.16 

19.7!oi 
20.81^ 


9(6) 
8(6) 
8(6) 


15.4(8.4) 
18.7(13.4) 
21.5(15.8) 



Table 1 resumes the results of our analysis in the 3 P-channels. In particular values of n 
and Q 2 outside (inside) parentheses refer to the option N — 1 (2) in eq. ([|). The values of 777-2 
are above the threshold of open charm (3727 MeV), therefore according to our predictions 
they are resonant states, which decay into charmed particles. We observe that the above 
mentioned self-consistency condition is satisfied for all 3 P-channels. Another condition, 
necessary in order to check stability of minima of Di with respect to parameters, consists 
in examining the two graphs corresponding to eqs. mi(m°, s ) = mcs, 777 2 (7n°, s ) = mpR, 



where GS means "ground state" and FR "first resonant state". Such graphs, represented 
in figs. 4 and 5, are quite similar to those drawn for the vector channel [p}] , to the difference 



that in the vector channel both mas and mpa are known experimentally, whereas in the 
other channels tufr is deduced from our numerical procedure. The analogous behaviour 
witnesses the self-consistency of our method. This is why in the 3 P!-channel we have not 
chosen the lowest minimum of D\. Indeed this minimum - which yields erratic values of 
some parameters, like £j and sq, and a mass iri2 greater than the one of the first resonance 
of the 3 P 2 -channel - corresponds to a pair of graphs quite different from the others, as one 
can see in fig. 5b. 

Our mass predictions differ from others deduced on the basis of potential models [IB]; on 
the contrary they fulfill, within errors, the center-of-gravity rule derived for P-channels of 
Charmonium (see|§ and refs. therein). 

We have already shown [p] , |l~2|j that g\ contains useful information for predicting some 
partial decay widths. The formulae for decay withs of Xco(2) into two gammas (gluons) read 

^ " 3(2J+l)m UJ ' X ^ 9 " 8 UJ r ^ 27 (7) 

where m = mcs, J (0 or 2) the spin of the resonance, 

9 2 = i^fjf> T = ^Tr[(f ! -4 + m e )0(} + m e )0), (8) 

k is the four-momentum of the quark, I an integer, O an operator and / an O-dependent 
normalization factor. For 3 Po(2) channels T turns out to depend critically on the modulus 
squared of the relative four- momentum p = 2k — qoi the quark with respect to the antiquark. 
The average of p 2 over the hadronic state results to be m 2 — Am 2 . According to Novikov 
et al. (see also ref. I12H ), we assume p 2 = —m 2 , therefore m 2 = |m 2 . In the 3 Po-channel 
/ Xc0 = 1, moreover O Xc0 is the identity operator, yielding / = 1 and T Xc0 ~ 0.268. As to the 
3 P2-channel, l Xc2 = 2 and 



°xc2 = ipPu + ivVn + , Vfivh = -Af - g^u, (9) 



2 / _<hQ. 

3 q 

Hence T Xc2 = and / = 1/5. Calculations yield 



r 

L XcO- 


^2 7 = 


14.lt^ 4 (5.6^ : 4 ) keV, 


(10) 


r 


+19 = 


26.81?:! ( 14 -° ± 5 ) MeV , 


(11) 


r 

1 Xc2~ 


^2 7 = 


0.99 ± 0.09 (0.32±°;^|) keV, 


(12) 


r 

L Xc2~ 


-2g = 


1.88±g;ii (1-73 ±0.16) MeV, 


(13) 



parentheses reporting experimental values ]p"7f|. 

As regards Xci, we assume, according to Novikov et al.||, the decay into light hadrons 
to be dominated by the two-step process xa ~^ 9 9*, 9* - ► QQ, i- e -> 



= 7^7 C ^™ s ? / d<bz- -h, u H^, (14) 



1 9 4, x „ r , 11 
2m x 3 J g 4 4 

where <i$ 3 is the phase space of the gqq system, h^ u (H^) the hadronic tensor of light 
(charmed) quarks and q 2 the effective mass squared of the qq system. Furthermore ~g Xcl is 
given by eqs. (§), with l Xcl = 2 and Xcl = 757^, whence T Xcl = 1 and / = 1. Then eq. (|i~4T) 
yields 

4 12m 6 a 3 , „ f - 5 , 2 1 - cosOcoscj) 



Ty 1^000 = t: — o — 5-^ / cfcos# / dx UJ x —. , (15) 

where e = 5/[l + 5 — cos9(l — 5)}, 5 = ql/m 2 , q 2 is a lower limit to q 2 , u> the energy of 
the "real" gluon, 9 (0) the angle between the direction of the real gluon and the one of the 
(anti-) quark in the overall cms. q 2 , uj and <fi depend on 9 and x through four- momentum 
conservation constraints. Taking into account the value of g\ (see table 1), and assuming, 
according to the lower bound of perturbative QCD scales, q 2 = 4 - 5 GeV 2 , we get Y Xcl ^ gq q 
= (658^207) keV, quite consistent with the experimental data, (640 ± 100) keV. 



Other authors |Tq, 19, 20] calculate partial decay widths of Charmonium P-wave reso- 



nances by means of the factorization theorem, including colour octet contributions in order 



to eliminate infrared divergences Ell]. They do not calculate the nonperturbative factors 



involved, rather they determine them either from some of the partial decay widths or from 



low-energy e + — e~ data|22|]. On the contrary, we do not take into account the octet contri- 
bution (relatively small in percentage); however, thanks to QCD sum rules, we can calculate 
nonperturbative factors. 
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Figure 1: 3 P channel: behaviour of m x vs n for £ = 2.5 
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Figure 3: 3 P 2 channel: D 2 vs g±. The arrow indicates the oblique inflexion. 
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Figure 4: 3 P and 3 P 2 channels: m° c vs s . The full lines refer to eq. mi = m GS , the dashed 
ones to eq. m 2 = m FR) where GS(FR) means "ground state" ("first resonant state"). 




Figure 5: 3 Pi channel: mj! vs Sq. See caption of fig. 4. Graphs b) refer to the last minimum 
of D\ and present a clear anomaly with respect to the others. 



